! IMPORTANT NOTE
! These weight need to be multiplied by 1/2
        hx(dcellx  ) = 0.5_num + cell_frac_x
        hx(dcellx+1) = 0.5_num - cell_frac_x
